astro-ph/OOxxxxx, KNGU-INFO-PH-15 
Enhanced black hole formation by multiple winding cosmic strings 

Michiyasu Nagasawa 

Department of Information Science, Faculty of Science, Kanagawa University, Kanagawa 

259-1293, Japan 
Electronic address: nagasawa@info.kanagawa-u.ac.jp 

Abstract 



Primordial black hole formation by cosmic string collapse is reconsidered 
when the winding number of the string is larger than the unity. The line 
energy density of the multiple winding string becomes greater than that of 
the single winding string so that the probability of the black hole formation 
by the string collapse during the loop oscillation would be strongly enhanced. 
Although the production of the multiple winding defect is suppressed and its 
number density should be small, the enhancement of the black hole formation 
by the increased energy density may provide a large number of evaporating 
black holes in the present universe which gives more stringent constraint on 
the string model compared to the ordinary string scenario. 
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I. INTRODUCTION 



Topological defects [|I| are considered to be produced in the early universe during the 
cosmological phase transitions accompanied by some kinds of symmetry breaking. They 
would contribute not only to the experimental investigation of particle physics models in 
our universe as a high energy laboratory but also to the theoretical explanation of various 
unresolved problems in the standard Big-Bang cosmology. Here we consider the primordial 
black hole formation by cosmic strings which is one of the notable subjects both in the 
observational astrophysics and the modern particle cosmology. 

A cosmic string is a line-like topological defect which can be characterized almost by one 
model parameter, 77, the symmetry breaking scale of the phase transition in which the string 
is generated. Particularly the line energy density of the string, /i, is expressed as 

/i ~ , (1) 

using 77. In the course of the cosmic evolution, the size and number distribution of the string 
is believed to obey the scaling model. In the scaling distribution, the ratio of the energy 
density of the cosmic string to the cosmological background energy density should be time 
independent. The numerical value of this ratio can be written approximately by 

where Mpi is the Planck mass. Note that when rj ~ the grand unification scale, that is, 
lO^s-i^ GeV, 

G/i ~ 10"^ , (3) 

which is suitable for the initial seed amplitude of the cosmological structure formation. This 
is why the GUT scale string has been not only regarded as one of the possible consequence 
from some unified theories but also considered to be a promising candidate of the required 
density perturbation seed for galaxies and clusters of galaxies. The viability of such a 
scenario can be investigated by the direct comparison between the results of numerical 
simulations and the various observational facts such as the galaxy distribution and the 
cosmic microwave background radiation fluctuation similarly to other cosmological structure 
formation scenarios. Although in the literature, the string scenario seems to be a less 
promising one than that based on the inflationary epoch, the fact that there is no established 
version of cosmological structure formation scenario suggests cosmological defects may have 
played an important role in the cosmic evolution history. 

In this paper, we pay particular attention to the constraint on the model parameter in 
which the cosmic string is produced using the black hole formation by loop oscillation. The 
point which has not been discussed so far is we employ the multiple winding number string as 
black hole sources. In the next section and the third section, we briefly review the primordial 
black hole formation by cosmic strings and the production of multiple winding strings. Based 
on the results obtained in these two sections, the revised observational constraint on the 
string model is shown in the section 4. The final section is devoted to discussions and 
conclusions. 
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II. PBHS FROM COSMIC STRINGS 



In this section, let we briefly summarize the primordial black hole formation from os- 
cillating string loops. Black holes are one of the most exciting consequence of the general 
relativistic gravitation theory and their astrophysical and cosmological influences have been 
investigated by many people. Although the conventional source of their production is the 
gravitational collapse of astronomical bodies, the primordial formation in the early universe 
is also an interesting and important issue which is worth to be studied in detail. 

Particularly black holes whose mass is M,, ~ lO^^g have the life-time which is comparable 
to the age of the universe so that they would evaporate at present and may provide sources 
of ultra high energy cosmic rays including extragalactic 7-rays [0 and other astroparticle 
phenomena whose mechanisms have not been solved completely. In addition to that, the 
relic black holes might contribute to the energy density of the present universe and solve the 
dark matter problem which is one of the most important cosmological problems . Such 
primordial black hole formation could be predicted not only by the density perturbations 
produced during the inflationary expansion of the universe but also by the collapses of 
oscillating string loops. 

In the standard picture of the cosmic string evolution, it is believed that the distribution 
of the size and number of strings obeys the one scale model, that is, the scaling distribution. 
In this formula, the string can be classified to two kinds of string, one is an infinitely long 
horizon-scale string and the other is a loop. As we mentioned in the previous section, the 
relative energy density of infinite strings is scale-invariant since they lose their energy in 
the course of cosmic evolution. In the conventional model, the main channel of the infinite 
string energy loss is considered to be a loop generation by the small-scale structure of long 
strings. Although there exists a recent claim which says that the infinitely long strings have 
very small wiggles on them so that the dominant energy loss mechanism of the string should 
be a direct particle production Q], the conclusion that the resulting string distribution can 
be described by the scaling formula is unchanged. In this formula, the energy density of the 
infinite string can be written as 

Poo = , (4) 

where L is the characteristic length scale of the string network which can be defined as 

L={CVt, (5) 

using the parameter, C, which is a typical number of infinite strings within the horizon. On 
the other hand, the loop number density n(/,t) for the loop size, /, and the time, t, can be 
expressed as 

where v is related to the loop creation length, 7^, and C, as 

^^CVl ^ (7) 
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and r is the efficiency of the gravitational wave radiation. If other way of energy emission 
is dominant, FG/x should be replaced appropriately. 

Whether the energy loss mechanism of the string loop is the gravitational wave radiation 
or not, oscillating loops may collapse to black holes. When the loop whose total length is 
equal to I is compactified into the region whose radius is 2Gfil during its oscillation, all the 
energy the loop has is contained within the gravitational radius so that a black hole of mass 
nl would appear. Although such a procedure still remains a conjecture, it is believed that 
the probability of black hole formation during a loop oscillation, P^h^ can be written as 

Phn ~ {Gf^r , (8) 
where k is a numerical coefficient and a is a power index which have been estimated as 

a ~ ~ 4 , (9) 

or [I 

a ~ -0.5 ~ 0.5 , (10) 

and 

a = 4.1 ±0.1, (11) 

n = 10^-9±°-2 , (12) 

based on various assumptions and calculations. Hereafter for simplicity, we employ a nu- 
merical value a = 4 as the most distinct case. Since k is cancelled out in the calculation 
process, its value does not affect the final conclusion. 

Combined with the observational constraints obtained from the existence of primor- 
dial black holes, the upper limit on the important string parameter, Gfi, can be derived. 
Although the evolution and the final destiny of the black hole have not been completely 
determined yet, there are two kinds of constraints by the evaporation of the black hole and 
its massive relics, which are naturally deduced from the assumption that the black hole 
loses its energy and mass due to the Hawking radiation. The former constraint is obtained 
from the observation of galactic and extra-galactic cosmic rays such as 7-ray bursts, 7-ray 
background, anti-protons and so on. The latter is based on the consideration that the black 
hole should not disappear away but leave a massive relic. Then it would contribute to the 
energy density of the universe and the total mass can be limited by the fact that the cosmic 
evolution must obey the standard history and the universe must not be over-closed. 

Although there are many uncertainties such as what the final state of the black hole is, 
how much the loop number density is and how much the black hole formation probability 
is, the summarized constraint results in ||^ 

G/i < (1 - 3) X 10"^ . (13) 

Note that this is in rough agreement with the bound by the cosmic microwave background 
radiation anisotropy as 

G/i < (1 - 2) X 10"^ , (14) 

and the most stringent constraint on the string model by the pulsar timing observation ||10| . 
As we mentioned in the previous section, the above conditions become marginal when rj ~ 
the GUT scale. 
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III. MULTIPLE WINDING STRING 



In the literature, only the simplest string configuration, that is, a string whose winding 
number, n, equals one was investigated since it was implicitly assumed that it seemed to 
be natural. Recently, however, it is claimed that the multiple winding string may play an 
important role in the cosmological evolution of the early universe. One is the initial condition 
for the infiationary universe and the other is the electroweak baryogenesis. 

The infiation scenario is the most promising paradigm which can solve many problems 
the standard Big-Bang theory cannot explain. So far, however, there has been no convincing 
model of the infiation and people have invented many models. Among them, the topological 



infiation |TT[ model can naturally provide an initial condition for the infiationary expansion 
since at the core of the topological defect the infiaton energy takes large amplitude which 
enables the dominance of the effectively constant energy density. The difficulty of this 
efficient model is that in order to satisfy the condition that the whole universe within the 
horizon must be dominated by the infiaton field, the symmetry breaking scale of the defect 
forming phase transition should be sufficiently large as 

V^Mpi. (15) 

Since no one knows a reliable physical theory at such a high energy scale, it would be useful 
if this constraint can be lowered. In the case of the topological string, the multiple winding 
number can settle the situation since the larger the winding number becomes, the larger the 
core size of the string grows, which results that the required value of t] decreases. Particularly 
when 

n = 3 X 10^ , (16) 



the topological infiation can occur at the GUT scale |12 |. 

The baryon asymmetry problem is one of the most significant problems in modern cos- 
mology and particle physics. Observations show the fact that there exists small baryon 
number as 



^ ~ 10-1° , (17) 



which is the ratio of the baryon number density to the entropy. To explain the generation of 
baryon asymmetry in the course of cosmic evolution, three necessary conditions are required. 
These are the baryon number violation, C and CP violation and the deviation from the 
thermal equilibrium. The first and the second conditions may be satisfied by the appropriate 
model construction of particle interactions. On the contrary, the last one should be provided 
during the dynamical evolution of the universe. 

Because of the sphaleron transition process, all the baryon asymmetry should be erased 
at the electroweak epoch unless there exists an asymmetry between baryon number and lep- 
ton number. This is the reason why the electroweak baryogenesis has been regarded as the 



most promising scenario of baryon number generation for many years |13]. In the conven- 
tional model of the electroweak baryogenesis, the deviation from the thermal equilibrium is 
achieved by the propagation of nucleated bubbles and the interaction between bubble walls 
and the surrounding plasma during the first order electroweak phase transition. However, in 
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the standard model and its simple extension, the degree of the electroweak phase transition 
seems to be too weak to realize the first order transition. 

In such a situation, various alternatives have been suggested and the electroweak string 
has been proved to bring the non-equilibrium condition even if the electroweak phase tran- 
sition is not of first order since it should be left after the transition and its collapse resem- 



bles the bubble wall propagation |jT^. Although the idea that the defect can realize the 
out-of-equilibrium remains valid, it has been shown that the initial production rate of the 
electroweak string in the standard model is too low so that it is difficult to explain the ob- 
servational value quantitatively by this scenario [|r^ . Recently more effective scenario using 
the string has been proposed ||T6|. In this case, when 

n > 2 ~ 3 , (18) 

the deviation from the thermal equilibrium can be satisfied by sphaleron bound states on 
strings and their following decay. 

Therefore it would be useful if the existence of the multiple winding string can be con- 
strained by some astronomical and cosmological observations and here we consider the pri- 
mordial black holes produced by the string loop oscillation. In order to perform a quantita- 
tive analysis, we estimate the formation probability of the multiple winding string. Although 
similar consideration can be applied to the multiple winding monopole which may be useful 
for the topological infiation, hereafter we will concentrate on the string case. 

We employ the Abelian Higgs model with the Lagrangian as 

^ = -If'^'^F,^ - 1{D^<I>)\D,<P) - iA(0t0 - , (19) 
D^ = d^- leA^ , (20) 
F^u = d^A, - d,A^ , (21) 

where is a complex scalar field, A^ is a gauge vector field, e is the gauge coupling constant. 



In this model, there is a well-known string solution called the Nielsen-Olesen vortex |T7 
which is a two-dimensional slice of the string and the string configuration with a fixed 
winding number, n, can be determined by one parameter, /?, which can be defined as 

/3 = 4=(-r. (22) 

where iris is the mass of the scalar field and rriy is the mass of the vector field. The important 
characteristic of the string configuration is that as j3 or the winding number, n, increases, 
the width scale of the string core also increases. 

The stability of the multiple winding string has been analyzed and in some parameter 
range, it is stable. For example, the calculation of the string line energy density shows 
that the multiple winding string is stable for the case (3 < 1 and unstable for (3 > 1 [jl8| 



Moreover, the interaction force between two vortices is attractive when /3 < 1 and repulsive 
when /5 > 1 0]. 

Now we estimate the formation probability of the multiple winding string. In this paper, 
we review the calculation process briefiy and the detailed method can be found in the 
reference ||20|| . 
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Before introducing how to estimate the string formation probabihty for the multiple 
winding, let us summarize that for the single winding case In the Kibble mechanism, 
the phase of the Higgs field is considered to take random value at each correlated region 
whose size is of scale, ^, which is the correlation length of the fields. The result shows that 
the formation probability of the single winding string, P(l) is 

^(1) = ^/2vr = \ . (23) 

In this procedure, the string identification process can be described as follows. First divide 
the physical space into domains whose size is typically of ^, that is, the correlation scale. 
Then assign the phase of the Higgs field randomly to each representative point of every 
domain. Finally interpolate the phase between two representative points so that the gradient 
energy of the Higgs field takes the minimum value. Thus we can distinguish the region where 
there is a winding number and a string exists from that where there is no winding number 
and we cannot find a string. 

Based on the above argument let us proceed to the multiple winding string version. If 
we simply apply the Kibble mechanism procedure to the multiple winding case, we have 
to encounter a problem. Since the phase difference between two neighboring representative 
points in the triangle division method which we have described, the total difference of the 
phase along the triangle circumference should be smaller than or equal to Svr, which means 
it must be 27r at most because of the phase continuity. Thus the winding number cannot be 
larger than the unity. The most natural solution to this situation may be the employment of 
other polygon than the triangle, for example, a hexagon. In the hexagon case, the total phase 
difference is improved to be Gvr and the region within which the winding number becomes 
two can be allowed. However, there remains a problem even in this case because we cannot 
distinguish the state there is one piece of double winding string from the state there are two 
pieces of single winding string. Moreover, if we apply the triangle method to the area where 
one double winding string can be found by the hexagon method, we will misunderstand that 
there is one single winding string. Therefore not only the winding number but also even the 
magnitude of the total winding depends on the method of space dividing. 

Thus in order to estimate the formation probability of the multiple winding string cor- 
rectly, we introduce a new scale, i?, within which there should be only one piece of string. 
Then the number of vertices the polygon has becomes ttR/^. Thus the total phase difference 
is 

ttR tt^R 

The most natural interpretation would be that R is comparable to the diameter of the string 
core which is approximately equal to the inverse of the scalar mass which results 

^ = /rJ rr N ~ ' (25) 

ms{T = Tg) XTg 

where Tg is the temperature at which the phase transition terminates and the configuration 
of topological defects cannot be erased by the thermal fluctuations of fields. On the other 



hand, the correlation scale, ^, can be expressed as 
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(26) 



As a consequence, the possible maximum winding number which belongs to one string should 
be 







TT 


2n^_ 




2A 



(27) 



where [ ] denotes the Gauss's symbol. Thus we can say that the multiple winding string can 
be produced when A << 1. 

Even in the case A ~ 1, we cannot completely deny the formation of the multiple winding 
configuration since in the actual situation, the correlation length scale can be inhomogeneous, 
that is, be various at each correlated region. Based on the result of the numerical calculation 
using the toy model [Q, we estimate the black hole formation possibility in the next section. 

Before closing this section, let us mention another possibility. When the self-coupling 
constant of the scalar field is smaller than the gauge coupling constant, that is, /5 < 1 for the 
local gauged string case, the attractive force acts on strings and they would accumulate to a 
multiple winding string. In this case, the formation probability of the string whose winding 
number equals n, P{n), can be easily calculated as 



Pin) 



P(l) 



n 



(28) 



This is also the case which the black hole formation probability is calculated in the following 
section. 



IV. PBHS FROM MULTIPLE WINDING LOOPS 

In this section, we estimate the primordial black hole formation probability by the os- 
cillating string loop and the resulting modified observational constraint. Since we have the 
formula for the ordinary single winding string case in (|13D, the only thing we have to do 
is to clarify the difference when we employ the multiple winding one. The constraint (|T3p 
was deduced by the upper bound on Qbh, the relative energy density of black holes to the 
critical energy density of the present universe. Hence in the following we will derive that for 
the n winding string, Q^^- 

Here we assume the modified formula of the parameter which characterizes the multiple 
winding string. If the line energy density of n winding string, be written as 

/i„ = A{n)iJ. , (29) 

using the enhancement factor, A{n) > 1, then 

A{n) < n , (30) 

since otherwise a multiple winding string becomes energetically unstable to the division into 
n pieces of single winding strings. As we have introduced in the equation the formation 



8 



probability of the black hole by the loop oscillation is in proportion to the power of the 
string line energy density so that this probability can be enhanced as 



^BH 



A{n) 
Ml) 



A{nr ^ h , (31) 



where the superscript of Pbh depicts the corresponding string winding number and we define 
this enhancement factor as /i. 

In order to see how the constraint on the string model should be replaced, it is also 
necessary to estimate the seed loop number density in addition to the above black hole 
formation probability. Using the assumption in the equation (^9]) and one more assumption 
that the multiple winding string also obeys the scaling distribution described by the formula 
(^), we can calculate explicitly two modifications by the increase of the winding number 
which affect the black hole number density at present. One is the reduction factor which 
comes simply from the formation probability of multiple winding strings and the other is 
the enhancement factor which is related to earlier formation of black holes. 

First the energy density of n winding strings, p^, is smaller than that of single winding 
strings, Poo, as 

pl. = P{n)^,. (32) 

where we have written the probability of n winding string production per correlation volume 
at the initial formation epoch as P{n) similarly to the previous section. Then there appears 
a modification to the final formula of the black hole density parameter as 

oc P{n) , (33) 

which results an relative difference compared to a single winding 

f - ^^^^ r34i 

h = ^y (34) 

where we define this reduction factor as f2- 

The next point is the time, f^ui ^^^t the black hole of specific mass M^, which is evapo- 
rating at present was formed becomes earlier in the multiple winding string case. Since 
can be calculated by the relation as 

= fin^tsH , (35) 

t^fj becomes smaller as /i„ increases. Strictly speaking the value of M* should depend on the 
time of formation. However, the rate of change of the black hole mass, M, is proportional 



to the power of M as p3 



^ « . (36) 

Then the life-time of the black hole is much larger than its creation epoch. For example, 
~ 10~^^ sec for the GUT scale string is sufficiently smaller than the age of the universe. 



9 



to ~ 10^^ sec. Thus it is a good approximation that is independent of the string winding 
number, n, so that we can say that 

oc - . (37) 



One more note is that in the equation ( p5|) it is assumed that the black hole is formed 
instantaneously after the loop creation. However, since the important feature is the relation 
between t^^^ and we do not investigate whether this is appropriate or not in detail. 

In order to calculate fi^/^, we have to take various mass of the black hole into account. 
Hereafter, however, we only estimate the number density of the black hole whose mass is M^,. 
This simple consideration is sufficient to see the essential dependence on /i„ of ^"bh since 
the black hole of mass M^, has the most dominant effect and the exact treatment produces 
the completely same result. The initial black hole number density whose mass equals M*, 
n%jj{t'%fj), can be calculated as 

= Pbh'^7oop{'^^h) ■> (38) 
where n^oop is the loop number density of the n winding string and can be written as 



nlopitln) - ^^^^ , (39) 



where Ap^ is the energy loss of an infinitely long n winding string to loops. Using the 
formula (^), 

^Ploitln) ^ Cp-^{tl^) = CP{n)C-^, , (40) 

where C is the efficiency of the loop production. When we leave only the dependence on 
the winding number, 

uIMh) ^ PBHPin)-^ . (41) 

['^bh) 

The following evolution of the black hole number density, n'^j^, obeys the ordinary dust 
particle law as 

r^S^(t)(xa(t)-^ (42) 

where a{t) is the cosmological scale factor. Then since t^^j is in the radiation dominated 
era for the parameter region we are interested in, the present number density of the black 
hole whose mass is equal to can be expressed as 



a{t 



BHj 



a(t( 

teq J V to 



j-n \ 3/2 , 2 



nlnitBH) C-^] , (43) 
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in the flat universe where teq is the radiation matter equahty time. Then the density pa- 
rameter for the black hole can be written as 



Po 



(44) 



if we neglect the mass loss of the black hole where po is the critical density of the universe. 
Combined with the relations (^) and (pTj), we can obtain the flnal relation as 

^BH P'^HP{n)iiT , (45) 



and we deflne an additional modiflcation factor, f^, as 



,3/2 



(46) 



As a result, the overall modification factor of string whose winding number 

is equal to n, f{n), should be described as 

fin) = hhh = ^^Ainr^/^ , (47) 

using the equations (|3TD, (|3^) and (^61). Thus when A{n) > 1 at a certain n, this enhance- 
ment may overcome the reduction by the smallness of P{n) and make the present black hole 
number larger compared to the single winding string case although the loop creation size 
and other parameters are not completely determined in the multiple winding string case. 

At last we can estimate the modification factor of the density parameter of the black hole 
at present for the multiple winding string. The first case is the string formation probability 
is calculated including the fiuctuation of the correlated region size. The assumption for the 
parameter is 

a = 4 , A{n) = n , (48) 

which makes the modification rather large. The final expression of the factor, f{n), can be 
written as 

P(n) 



fin) 



Pil] 



-n 



5.5 



(49) 



and the numerical values for n = 1 — 4 are depicted in Table 1. In this case, /i and /s is too 
small to enhance the black hole production so that the result is trivial, that is, the constraint 
on the string model parameter is unchanged since the most dominant contribution to the 
7— ray emission and other astrophysical phenomena would be provided by single winding 
strings. Although we do not know the numerical value of P{n) for much larger n, it would 
be natural that f{n) is the decreasing function of n even in such cases. 



n 


P{n) 


fin) 


1 


0.2102 


1 


2 


8.36 X 10-^ 


1.80 X 10-1 


3 


4.8 X 10-^ 


9.6 X 10"^ 


4 


8 X 10-11 


8 X 10-^ 



Table 1. Modification Factor for Multiple Winding String 
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In contrast to that, the other case in which the string coalescence is taken into account 
may be more interesting. Using the same assumption for the parameters, a and A(n)^ as 
the former case, the formula of the modification factor can be calculated as 



fin) 



n 



4.5 



(50) 



where we have substituted the equation (PSI). The factor f{n) are shown in Table 2 for 
n = 2 — 4. It can be obviously seen that the degree of the enhancement becomes greater 
as n increases. Therefore it may be plausible that the string which has a large number of 
winding might impose a more stringent constraint than the single winding string scenario. 



n 


2 


3 


4 


fin) 


22.6 


140 


512 



Table 2. Modification Factor for String Coalescence Case 



Finally we show the modified constraint formula for the string model parameter, t], which 
is significant from the point of view of the particle cosmology. As we have seen in this section, 
the symmetry breaking parameter is related to the black hole number as 



a+3/2 



(51) 



Then the revised upper bound on Gt] for the multiple winding string case, Gr] can be 
expressed as 



fin) 



-2/(3+2a) 



Gt]"- 



(52) 



where Grj'^l^ is the upper bound on Gr]"^ for the case in which the multiple winding string 
is not incorporated. 



V. DISCUSSION AND CONCLUSION 

In this paper, we have estimated how the observational constraint on the particle physics 
model parameter in which the cosmic string is produced should be modified when we con- 
sider the multiple winding string which may be useful for the topological inflation and the 
electroweak baryogenesis. The modified formula of the constraint on the line energy density 
which can be translated to that on the symmetry breaking scale, rj, can be written as 

< /(n)-2/(3+2«) X 10-6 , (53) 

in which we assume the simple scaling distribution of the string. 

The numerical values of the above modification factor are calculated in two cases. One 
is the case that the multiple winding string is produced at the initial formation epoch and 
the modification is not so significant since 

fin) < 1 . (54) 
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In the other case, the dynamical evolution of the string after its formation is concerned and 
the string accumulation by the attractive interaction between strings makes the enhancement 
remarkable as 

f{n) » 1 , (55) 

which means the upper bound on rj must be much lower. 

Thus we can say that a more stringent constraint by primordial black holes on the model 
parameter might be obtained by multiple winding strings compared to the conventional 
scenario of ordinary single winding strings. Just for the reference, let us estimate the case 
in which all the winding number within the horizon are concentrated on a single string. At 
the string formation epoch, the number of the correlated region in the horizon volume can 
be estimated as 

^=(^)'- (56) 

When we take the annihilation of a string by an anti-string into consideration, the bulk total 
winding within the horizon would be 



'P{1)N . (57) 
If we assume T rj, then 

using the relation between t and T in the standard cosmology where Af is the gauge degree 
of freedom of particles at T. Then the modified constraint on rj can be written as 

/ X l/2+6/(2a+3) / /T? \ 3(2Q-l)/2(2a+3) 

'"^ 4^ 10-'. (59) 



yM,,j ~ \a.3P(,iy 

leaving the parameter a in the formula (^. When a = 4, the resulting constraint can be 
expressed as 



V 



M 



pi 



< 10-^ , (60) 



which is marginally compatible with the GUT string, rj ~ 10^^ GeV and a little more rigorous 
than the original one (^3|). It might be the most useful condition in order to constrain the 
model construction of the string generation in the early universe even if the cosmic microwave 
background radiation fluctuation and the pulsar timing are included. 

This is not the end of the story. There may be other effects which can enable the 
abundant formation of the multiple winding defect such as the relaxation of the geodesic rule 



2J] which is applied when the phase between two different correlated regions is interpolated. 
Moreover, it may be possible that the homogeneous gauge flux distribution greatly enhances 
the total winding number within the correlated domain so that huge winding strings might 



be produced by the string coalescence |25|. Further work would be needed in the future. 
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